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D . Abstract 

Ti^lj- ! We compute the magnetic dipole transitions between low lying Heavy Quarkonium states in a 



Ph. 



model independent way. We use the weak-coupling version of the effective field theory named 
potential NRQCD with the static potential exactly incorporated in the leading order Hamiltonian. 



jil^i The precision we reach is k^/m? x 0{a1,v'^) and k^/m? x 0{v^) for the allowed and forbidden 

1-^1 transitions respectively. We also resum the large logarithms associated to the heavy quark mass 

scale. The specific transitions considered in this paper are the following: T(15) — t- 775(15)7, 
. J/i^{lS) ^ r?c(15)7, h,,{lP) ^ xw,i(l^)7, Xh2{lP) ^ hf^ilP)^, T{2S) ^ %(25)7, T{2S) 

10 , rH){lS)^ and rjh{2S) — )■ T(1S')7. The effect of the new power counting is found to be large and 

rn 

' the exact treatment of the soft logarithms of the static potential makes the factorization scale 

O ' 

ff^ ■ dependence much smaller. The convergence for the bb ground state is quite good, and also quite 

reasonable for the cc ground state and the bb IP state. For all of them we give solid predictions. For 

^ I the 25 decays the situation is less conclusive, yet our results are perfectly consistent with existing 

H ■ 

data, as the previous disagreement with experiment for the T(25) — t- %(15)7 decay fades away. 
We also profit to compute some expectation values like the electromagnetic radius, (r^), or (p^). 
We find (r^) to be nicely convergent in all cases, whereas the convergence of (p^) is typically worse. 

PACS numbers: 12.38.-t,12.39.Hg,13.20.Gd,12.38.Cy 
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I. INTRODUCTION 



The heavy quarkonium has always been thought to be the "Hydrogen atom" of QCD. The 
reason is that the heavy quarks in the bound state move at non-relativistic velocities: f <^ 1. 
This allows testing the dynamics associated to the gluonic and light-quark degrees of freedom 
in a kinematic regime otherwise unreachable with only light degrees of freedom. Effective 



field theories (EFT's) directly derived from QCD like NRQCD [1] or pNRQCD [2] (for some 
reviews see Refs. js, 4|) disentangle the dynamics of the heavy quarks from the dynamics 
of the light degrees of freedom efficiently and in a model independent way. They profit 
from the fact that the dynamics of the bound state system is characterized by, at least, three 
widely separated scales: hard (the mass m of the heavy quarks), soft (the relative momentum 
IpI ~ mv ^ m of the heavy- quark-ant iquark pair in the center of mass frame), and ultrasoft 
(the typical kinetic energy E ~ mv'^ of the heavy quark in the bound state system). 

In this paper we use pNRQCD. This EFT takes full advantage of the hierarchy of scales 
that appear in the system, 

m ^ mv ^ mv^ ■ ■ • , (1) 

and makes a systematic and natural connection between quantum field theory and the 
Schrodinger equation. Schematically the EFT takes the form 

ido- — -V^^\r)\<P{f) = Q 
m J 

+ corrections to the potential ' pNRQCD 

-|- interactions with other low-energy degrees of freedom 

where Vs^\r) is the static potential and </)(r*) is the Q-Q wave function. 

The specific construction details of pNRQCD are slightly different depending upon the 
relative size between the soft and the Aqcd scale. Two main situations are distinguished, 
named the weak {2!, Q [mv ^ Aqcd) and the strong [6] {mv ~ Aqcd) coupling version of 
pNRQCD. One major difference between them is that in the former the potential can be 
computed in perturbation theory unlike in the latter. 

It is obvious that the weak-coupling version of pNRQCD is amenable for a theoretically 
much cleaner analysis. The functional dependence on the parameters of QCD {as and the 
heavy quark masses) is fully under control and directly derived from QCD. The observ- 
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ables can be computed in well defined expansion schemes with increasing accuracy, and 
nonperturbative effects are exponentially suppressed compared with the expansion in pow- 
ers of as- Therefore, observables that could be computed with the weak-coupling version 
of pNRQCD are of the greatest interest. They may produce stringent test of QCD in the 
weak-coupling regime (but yet with an all-order resummation of powers of as included) and, 
precision permitting, are ideal places on which to accurately determine some of the parame- 
ters of QCD. Nowadays there seems to be a growing consensus that the weak-coupling regime 
works properly for t-t production near threshold, the bottomonium ground state mass, and 
bottomonium sum rules. In reaching this conclusion is crucial the proper treatment of renor- 
malon effects, which leads to convergent series, and the resummation of large logarithms, 
which significantly diminish the factorization scale dependence of the observable. Neverthe- 
less, even in those cases, the situation is not optimal. For some observables, even if getting 
a convergent expansion, the corrections are large, or in the case of the bottomonium ground 
state hyperfine splitting a two-sigma level tension between experiment (see Ref. [7|) and 
theory [sl, |9[ exists. 

In order to improve the convergen ce p roperties of the theory, the perturbative expansion 
in pNRQCD was rearranged in Ref. 10|. In this new expansion scheme the static potential 

was exactly included in the leading order (LO) hamiltonian^^ ^The motivation behind this 

reorganization of the perturbative series is the observation [11] that, when comparing the 
static potential with lattice perturbation theory, one finds a nicely convergent sequence to 
the lattice data (at short distances). Yet, for low orders, the agreement is not good and 
the incorporation of corrections is compulsory to get a good agreement. This effect can be 
particularly important in observables that are more sensitive to the shape of the potential, 
and naturally leads to consider a double expansion in powers of v and as{m), where v has 
to do with the expectation value of the kinetic energy (or the static potential) in this new 
expansion scheme. 

In Ref. 10[ this new expansion scheme was applied to the computation of the heavy 
quarkonium inclusive electromagnetic decay ratios. It was observed an improvement of the 
convergence of the sequence for the top and bottom cases. It was particularly remarkable 
that the exact incorporation of the static potential allowed to obtain agreement between 
theory and experiment for the case of the Charmonium ground state. This leads towards 
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the second motivation of the present study. The possible apphcabihty of the weak-couphng 
version of pNRQCD to the Charmonium (ground state) and the n = 2 excitation of the 
bottomonium. For those states the situation is more uncertain. Whereas in Refs. 
uM it is claimed that it is not possible to describe the bottomonium higher excitations in 

nnri 

perturbation theory, an opposite stand is taken in Refs. [9|, llSMlTj . We hope that we may 
shed some light on this issue as well. 

The above discussion basically refers to the determination of the heavy quarkonium mass 
and inclusive electromagnetic decay widths. Obviously there are more observables that can 
be considered. Some of those are the radiative transitions: H{n) — > H{n')'y, where n, n' 
stands for the principal quantum numbers of the Heavy Quarkonium. In Ref. [ISi] the allowed 
{n = n') and hindered (n 7^ n') magnetic dipole (Ml) transitions between low lying Heavy 
Quarkonium states were studied with pNRQCD in the strict weak-coupling limit. They also 



performed a detailed comparison of the EFT and potential model (see Refs. [19|, |20| for 
some reviews) results. The specific transitions considered in that paper were the following: 
J/^(15) ^ r7,(15)7, T(15) ^ VbilS)^, T(25) ^ Vbi^S)^, T(25) ^ 1/^(15)7, VbC^S) ^ 
T(1S')7, hb{lP) — > Xbo,i{^P)l and Xfe2(l-P) h^ilP)'^. Large errors were assigned to 
the pure ground state observables, specially for Charmonium, whereas disagreement with 
experimental bounds (at that time) was found for the hindered transition T{2S) — )■ 77b(lS') 7. 
In this paper we apply the new expansion scheme to those observables. The precision we 
reach is k^Jm? x 0{a1^v'^) and k'^Jm? x 0{v^) for the allowed and forbidden transitions 
respectively. Large hard logarithms (associated to the heavy quark mass) have also been 
resummed when appearing. The effect of the new power counting is found to be large and 
the exact treatment of the soft logarithms of the static potential makes the factorization 
scale dependence much smaller. The convergence for the hh ground state is quite good. This 
allows to give a solid prediction for the T(IS') — i- 77^(15') 7 transition with small errors. The 
convergence is also quite reasonable for the cc ground state and the hb IP state. For all of 
them we give solid predictions. For the Jjjp{lS) — ?■ ric{lS) 7 our central value is significantly 
different of the one obtained in Ref. |is| though perfectly compatible within errors. For 
the 25 decays the situation is less conclusive. Whereas for the T(25') — t- rib{2S)j decay we 
do not find convergence, previous disagreement with experiment for the hindered transition 
T(2S') — 7- 775(16') 7 fades away with the new expansion scheme. 
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The above observables depend on the expectation values of some quantum mechanical 
operators like, for instance, or f^^ (the electromagnetic radius). The study of them in 
an isolated way is interesting on its own. First, they provide us with a very nice check 
of the renormalon dominance picture. According to this picture the determination of the 
heavy quarkonium mass using the static potential (in the on-shell scheme) should yield a 
bad convergent series, as it is actually observed. The reason for this bad behavior is the 
existence of an r-independent constant that contributes to the potential and deteriorates 
the convergence of the perturbative series. If this is so, a check of this picture would be 
the computation of observables that are not affected by adding a constant to the potential. 
For those good convergence is expected. This is actually the case of (p^) or (r^). We 
nicely see in Sec. IIIII that this picture is confirmed. We find the electromagnetic radius 
(somewhat surprisingly) to be nicely convergent in all cases. This allows us to talk of the 
typical (electromagnetic) radius of the bound state in those cases. The kinetic energy is also 
(though typically less than the radius) convergent except for the 25* state. Then we can also 
define a typical velocity v = ^J {p'^) /m? for those states. 

The paper is distributed as follows. In Sec. [Tllwe discuss the theoretical background of 
the computation and display the formulas we use for the decays. In Sec. IIIII we analyze 
(p^) and (r*^), and discuss renormalon dominance. In Sec. IIVI we compute the radiative 
transitions. Finally, in Sec. |V]we summarize our main results and give our conclusions. 

II. THEORETICAL SETUP 

For the purposes of this paper we can skip most details of pNRQCD. We will only need the 
singlet static potential, V"i°^(r) — V{r), and the spin-dependent potential V'^^{r) — )■ V52(r)^. 
The static potential will be treated exactly by including it in the leading-order Hamiltonian 

H^""^ = -^ + V{r), and = ^n/0nz(r) , (2) 

where mr = 71111712/ {mi + 7712) (in this paper mi = m2 = m). The static potential will 
be approximated by a polynomial of order + 1 in powers of ctg (C/ = {N^ — l)/(2A'c), 

^ For simplicity we omit the index " s" for singlet and the upper index " (0)" and " (2)" throughout the paper. 
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Ca = N, 



N 



n=l 



477 



(3) 



In principle we would like to take as large as possible (though we also want to explore 
the dependence on A^). In practice we take the static potential at most up to N=3, i.e. 
up to 0{a'l) including also the leading ultrasoft corrections. This is the order to which the 
coefficients completely known: 



ai(z/; r) = Oi + 2[3q In {ue^'^r) 



IT 



a,{u; r) = a2 + + (4ai/3o + 2/3i) In (i/e-^^r) + 4^,' \n' (z/e^^r) , 

a^iu; r) = as + a.^^n^ + — /^oA + IQCsPo 

+ (^271^(3^ + 6a2/3o + 4ai/3i + 2/^2 + y^'ivr^^ In {ue^^r) 

+ (^12ai(3^ + 10/3o/3i^ In^ (z/e^^r) + 8(3^ In^ {ue^^r) 



(4) 



The C(q!s) term was 



term in Refs. 
piece in Refs. 



23 
26 



computed in Ref. 



2lj, the 0{al) in Ref. [2^, the 0{al) logarithmic 



, 24|, the light-flavour finite piece in Ref. 25|, and the pure gluonic finite 



27|. Explicit expression for can be found in Ref. {4]. For the ultrasoft 



corrections to the static potential we take 



(5) 



28 



We will not use the renormalization group improved ultrasoft expression in this paper 
as its numerical impact is small compared with other sources of error. 

We will always work with three light (massless) quarks. For the case of the bottomonium 
ground state we also incorporate the leading effect due to the charm mass: 



r = — - 



dx — I 1 



2x' 



-2mcrx 



(6) 



Its effect will be quite tiny. Therefore we have only incorporated Eq. (Q in our final (A^ = 3) 
evaluations and not considered any other subdominant effects in the charm mass. 
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The spin-dependent potential will be treated as a perturbation. It will contribute to the 
hindered Ml transitions. Nowadays this object is known with next-to-leading-log (NLL) 
accuracy 33|]. Nevertheless, for consistency with our accuracy, we will use its LL expression 

4 



(2) 



where 



3J] (see also 35j for the derivation in vNRQCD) 



depends on the NRQCD Wilson coefficients. With LL accuracy they read 

-Ca 



dsv{y) = dsv{m), 
4„(z/) = d^^{m) + 



Ca 



(3o - 20 A 



7ra,(m)(z'^o-2^^ - 1), 



where 



"o 1 / u 

~ 1 Qsi^) In I — 



dsvijn) = Cf ^C/ - 7ias{m 
dvv{m) = ■ 



(7) 



(9) 



(10) 



The theoretical study of the Ml transitions in the strict weak-coupling limit of pNRQCD 
has been carried out in detail in Ref. 



18| . A particular relevant result was that non- 



perturbative effects, associated to the mixing with the octet field, were subleading and 
beyond present precision. We can use their results in our power counting scheme with minor 
modifications (note that the dependence on the ultrasoft scale only enters marginally through 
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FIG. 1: Kinematics of the radiative transition H — t- H'^ in the rest frame of the initial- state 
quarkonium H . 



the static potential). The expressions we use for the decays are the following^ 

>2\ 



3 



_^ , 2\ 5 „./o(p^)nO 2 n'o(^g2(r))„o 

_5 / 2\ 5 „./o(p^)nO 2 n./o(Ks2(f))„o " 



r(n^Pj ^ n^Pi7) = — ^ — A — : — = -aej ^ 



2J + 1 

where in Eq. do = 1, di = 2, ^2 = 8/5 



{l + ny-d/-^ 



(11) 

2 

,(12) 

2 

, (13) 
(14) 



2M, 



H 



(15) 



and the anomalous magnetic moment of the heavy quark, which is renormalization group 



invariant, reads 



(16) 



In the following we use the notatfon {nS\p'^\nS) — (p^)no, {nP\p'^\nP) = (p^)ni, {n'S\p'^\nS) — „'o(p^)nO 
and so on. 
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K 



(1) = 



2n 



-(2) 



+ 



_31 57r2 
"l6 ^ I2" 



2 4 ' ^ 



(17) 
(18) 



317 



TT 



TT 



2ln2 



3C3 



144 8 



Ca 



2hnf 119 
36 ^ "36" 



TT 

"3 



We take k from Ref. [36] (it was originally computed in Refs. [37|, |38|, thought the first 
reference suffered from a factor 4 misprint). 



Eqs. flllfl2fl3|IT4|) follow from the expressions obtained in Ref. [18| except for the fol- 
lowing changes: (i) The matrix elements of r^, and V52 are computed using the exact 
solution of Eq. ([2]) with the static potential approximated to the power instead of using 
the Coulomb potential; (ii) we use the heavy quark anomalous dimension k to 0{a1); (iii) 
and our expression for V52, Eq. ([7]), incorporates the LL resummation of logarithms (this 
will actually be important for the 25* — )■ IS* decays). Overall our expressions are accurate 
with k^/m^ X (9(f^, a^) and k^/m? x 0{v'^) precision for the allowed and hindered transitions 
respectively, and also include the resummation of large (hard) logarithms. 

Eqs. ( fTT|12fl3lfT^ have been obtained in the on-shell scheme. Therefore they depend on 
the pole mass, m, and the static potential, V , both of which suffer from severe renormalon 
ambiguities. On the other hand, the decays themselves are renormalon free, as they are 
observables. Therefore, it is convenient to make the renormalon cancellation explicit. One 
first makes the substitution^ 



(m, V{r)) = {iTLx + 5mx-, Vx{r) - 2 5mx) 



where 



AT 



etc 



(19) 



(20) 



n=0 



represents a residual mass that encodes the pole mass renormalon contribution and X stands 
for the specific renormalon subtraction scheme. Matrix elements are renormalon free but 
not the heavy quark mass. Its renormalon ambiguity cancels with the one coming from 
the anomalous magnetic moment of the heavy quark. The renormalon structure of the 
chromomagnetic moment of the heavy quark has been studied in detail in Ref. \^'^. If one 



Note that 6mx and m (or V) have to be expanded to the same power in and at the same scale. 
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does the abelian-like limit one can get the renormalon structure of the anomalous magnetic 
moment of the heavy quark. One sees that it suffers from the very same renormalon than 
the heavy quark mass. Therefore the quantity (1 + K)/m is free of the renormalon ambiguity 
(or at least of the leading one). When rewritting the decay expressions from the on-shell 
scheme to the X scheme the change is absorbed in k so k — > where 

i^x = i^x^C(s{m) + K^^^al{m) + ■ ■ ■ (21) 



with 



(i)-^(i)_^5^(J)(^) (22) 



m " m 



^(2) _ ^(2) 



^5m«(^) - «:W^<5m(?)(^) + (^SmPi^))' . (23) 
mm mm \m m / 

Overall in Eqs. ( ITT1I12II13IIT41) we have to make the replacement (m, V, k) — )■ {mx,Vx, i^x) 

throughout. Note that once written in terms of renormalon free quantities one may consider 

different A^, A^' for vj^\ \ and the observable would still be renormalon free. 



III. APPLICABILITY OF WEAK-COUPLING TO HEAVY QUARKONIUM 

The allowed Ml radiative transitions depend on {p'^)ni and, at higher orders, on other 
expectation values such as {r'^)ni- To study them gives us a hint of the applicability of the 
weak-coupling version of pNRQCD to those states, and a very nice check of the renormalon 
dominance picture. In this section we compute the bound state energy (Eni), {p'^)ni and {r^)ni 
for the charmonium ground state and for n = 1,2 bottomonium states. In the cases where 
we have good convergence we will be able to obtain well defined values for Vni = a/ {p^)ni/'^'^ 
and {r'^)ni- 

Our reference values for the charm and bottom masses are mi,{mf,) = 4.19 \^ and 



fT^c{f^c) = 1-25 41 1 , w hich we then transform to renormalon subtracted schemes like the 
RS, RS' [42] or PS |43|. We will mainly use the RS' scheme and leave the RS and PS scheme 
for partial checking (in particular that the dependence on the renormalon subtraction scheme 
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is small). Therefore, we use (Sm^l, = and uj) = /3„/2^"*"^"' In 



Smil{^) = N^^S{l,b), 



5m 



(2) l^f^ 



6m 



(3) 



271 

/3o 



■5(M)^^+(|)5(2,6) 



=A^,., I — - I X 
/vr 



X 



where 



with Co = 1 and 



A r(n + l + &-fc) 
^(n,6) = X.^^ r(l + &-fc) 



+ 6- fc) 



2 ' 



Cl 



46/33 



-A 



and 



C2 



32/38 



(24) 



(25) 



(26) 



(27) 



6(6-1) 

For easy of reference, we give some typical values that we will use in this paper: 
"^6,RS'(0.7GeV) = 4902 MeV, mb,Rs/(lGeV) = 4859 MeV, me,RS'(0.7GeV) = 1648 MeV, 
'^c.RsKl GeV) = 1536 MeV. Our reference value for will be Nm = 0.574974 (for three 
light flavours) from Ref. 42|. To this number we will typically assign a 10% uncertainty. 
Our reference value for will be ai^^~^\l GeV) = 0.479778, which we obtain running 
down ai''^^^\Mz) = 0.118. We then run with four loop accuracy for the typical scales of 
the bound state system. Unless stated otherwise, throughout the paper we will set = i^f- 
The static potential we will consider in the following will be (in the RS' scheme) 

N 



n=0 

AT 

(V(^) + 26m'i'J)\,=yr ^ Yl ^RS',n«r'(l/r) if r < 



(28) 



n=0 



This expression encodes all the possible limits: 

(a). The case = 00, Uf = is nothing but the on-shell static potential at flxed order, 
i.e. Eq. Note that the = case reduces to a standard computation with a Coulomb 
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potential, for which we can compare with analytic results for the matrix elements. We will 
use this fact to check our numerical solutions of the Schroedinger equation. If we also switch 
off the hard logs our computation would be equal to the one performed in Ref. 18|]. We will 
use this fact to compare with their results throughout the paper. 

(b). The case Ur = oo (with finite non-zero Uf) is nothing but adding an r-independent 



constant to the static potential (see the discussion in Ref. pj|). Therefore, the result for 
{p'^)ni and {r'^)ni do not depend on the specific value of z/j (for a fixed heavy quark mass). 
In particular the value Uf = can be taken, which is equivalent to not considering any 
renormalon subtraction at all. On the other hand the binding energy Mni is renormalon 
dependent. This effect can be seen in full glory in Fig. [21 where we plot Miq and (r^)io for 
the case of the bottomonium using the static potential V^} at different orders in perturbation 
theory: = 0, 1, 2, 3. We clearly observe how for the i/^ = cxd, z/j = case, the bound state 
energy is not convergent (see dashed lines), whereas (r^)io is (see solid lines). On the other 
hand for the z/^ = oo, = 0.7 GeV case both the bound state energy and (r^)io show a nice 
convergent pattern as we increase (see solid lines). Note that (r^)io is exactly the same 
in both cases: z/j = or z/j = 0.7 (this is the reason only solid lines show up in Fig. Ob). 
The same analysis can be done for (p^)io, as one can see in Fig. [3] for the dashed lines (note 
though that (p^)io is less convergent than (r^)io). 

A rather similar picture is observed for the Charmonium ground state though, as expected, 
the sequences are less convergent. Again it is compulsory to incorporate the renormalon 
cancellation (finite Uf) to transform the bound state energy in a convergent sequence in A^, 
whereas (r^) and (p^) are always convergent (see the dashed lines of Fig. Hj). Everything in 
fully accordance with the renormalon dominance picture. 

For the n = 2 bottomonium states the situation is less conclusive. For both the P- and 
S'-wave (r^) is convergent, see the dashed lines of Figs. Oa and Oa respectively. On the 
other hand (p^) is only marginally convergent for P-wave (see the dashed lines of Fig. Ob), 
or even no convergent for the S wave (see the dashed lines of Fig. Ob), as each order is 
typically of the same size. 

(c). We now take z/^ = finite (and, for consistency, Ur > vj). We expect this case to 
improve over the previous results, as it incorporates the correct (logarithmically modulated) 
short distance behavior of the potential. Yet this has to be done with care in order not to 
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0(a) 



O(a^) 



O(a^) 



O(a^) 



0(a) 



O(a^) 



O(a^) 
O(a^) 



1.5 



2.5 3 3.5 
V (GeV) 



4.5 



FIG. 2: Plot of Mio and y^{r'^)io for the case of the hottomonium ground state using the static 
potential V^) at different orders in perturbation theory: = 0, 1, 2, 3. The dashed lines have been 
computed with Uf = 0. The continuous lines have been computed with Vf = 0.7 GeV. In both cases 
Ur = oo GeV. 
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1 L 



1 1.5 2 2.5 3 3.5 4 4.5 5 

V (GeV) 




V (GeV) 

FIG. 3: Plot of Yv^io o'^'^ ^lo /^^^ case o/ the bottomonium ground state using the static 
potential V-^^) at different orders in perturbation theory: N = 0,1,2,3. The dashed lines have been 
computed with = oo. The continuous lines have been computed with = 0.7 GeV. In both cases 
Uf = 0.7 GeV. 
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spoil the renormalon cancellation. For this it is compulsory from now on to keep a finite, 
non- vanishing, Uf, otherwise the renormalon cancellation is not achieved order by order in 
N, as it was discussed in detail in Ref. We have explored the effect of different values of 
Uf in our analysis. Large values of Uf imply a large infrared cutoff. This makes our scheme 
to become closer to a MS-like scheme. Such schemes still achieve renormalon cancellation, 
yet they jeopardize the power counting, as the residual mass does not count as mv'^. This 
has the price of making the consecutive terms of the perturbative series bigger. Therefore, 
we prefer values of Uf as low as possible, with the constraint that one should still obtain 
the renormalon cancellation, and that it is still possible to perform the expansion in powers 
of Og. In our analysis we observe that we can use a rather low value of Uf and yet obtain 
the renormalon cancellation. By also taking a low value of Ur we find that the convergence 
is accelerated and the scale dependence is significantly reduced. We illustrate this behavior 
in Figs. [31 m [5], and [6l where we can compare the case with Uj. = finite (continuous lines) 
and Ur = oo (dashed lines). This improvement is observed in all cases except for the 2S 
bottomonium (p^)2o- Specially relevant for us it that it accelerates the convergence of (p^)io 
for charmonium, and transforms (p^)2i into a convergent series. 

Leaving aside the 2S bottomonium state, it is particularly appealing to compare the = 
case for Ur = finite and z/^ = oo. The latter corresponds to the Coulomb approximation and 
it is the one used in the strict weak-coupling analysis performed for the radiative transitions 



in Ref. 



18(1 . One can see a very strong scale dependence, almost a vertical line compared 



with the Ur = 0.7 GeV case (see for instance Fig. [3]). Therefore, small variations of the 
scale produce very large changes in the theoretical prediction. This makes difficult to assign 
central values (and errors). This is not so after resumming the soft logarithms by setting 
Ur ^ 0. This produces fiatter plots. Note also that, typically, there is an scale where z/^ = cxd 
and i>r = 0.7 lines crosses. One can then take this scale as a way to fix the scale u of the 
computation with z/,, = oo (which correspond to the strict weak-coupling expansion). 

Overall, we find the electromagnetic radius (somewhat surprisingly) to be nicely conver- 
gent in all cases. This allows us to talk of the typical (electromagnetic) radius of these bound 
states. The kinetic energy is also (though typically less than the radius) convergent except 
for the 25" state. Then we can also define a v'^i = {v^) ni / 'fn^ for those states. We show 
these numbers in Tab. [H These numbers can be taken as estimates of the typical radius 
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FIG. 4: Plot of \J (r^)io and v\q for the case of the charmonium using the static potential at 
different orders in perturbation theory: N = 0,1,2,3. The dashed lines have been computed with 
Ur = oo- The continuous lines have been computed with z/,. = 0.7 GeV. In both cases Vf = 0.7 GeV. 
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FIG. 5: Plot of \J (r^)2i and w^i /or the case of the n = 2 P-wave bottomonium using the static 
potential at different orders in perturbation theory: N = 0,1, 2, 3. The dashed lines have been 

computed with v,. = oo. The continuous lines have been computed with = 0.7 GeV. In both cases 
uj = 0.7 GeV. 
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FIG. 6: Plot of \J (r^)2o o-nd v\q for the case of the n = 2 S-wave hottomonium using the static 
potential V^) at different orders in perturbation theory: = 0, 1, 2, 3. The dashed lines have been 
computed with v.,. = oo. The continuous lines have been computed with Uf. = 0.7 GeV. In both cases 
Vf =0.7 GeV. 
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of the bound state system and of the typical velocity of the heavy quarks inside the bound 
state. It is conforting that the numbers we obtain for f ^ are similar to those usually assigned 



either by potential models or by NRQCD (see for instance 4^, |45|). The specific values in 
the table have taken from the = 3 case at z/ = 1.5 GeV (and Uf = Ur = 0.7 GeV). For 
the bb ground state the result is very stable under scale variations, for the charm ground 
state and for the bottomonium P-wave the scale dependence is bigger. We stress that the 
numbers in the table should be taken as estimates. We do not attempt here to perform 
an specific error analysis of those numbers, as it is not needed for the decays. Let us just 
mention that one source of the error would come from the 1/m subleading potentials. In 
principle these effects would produce (9(f ^) corrections. For bottomonium and charmonium 
this would typically mean a ~ 7% and ~ 20% variations of the central values, respectively. 
Specially for bottomonium such uncertainties would compite with the difference between 
different N evaluations or, in some cases, with the scale variation. Finally we remark that 
for the 25* bottomonium state the numbers in the table should be taken with more caution 
as there is no convergence in the sequence in A^. One could actually be surprised by the 
fact that the (n, /) = (2, 1) and the (n, 0) = (2,0) states show this different behavior as far 
as convergence is concerned, since the typical transfer momentum is the same. One should 
note though that the (n, 0) = (2, 0) state has two nodes, and the most important one is a 
very low momentum. On the other hand this problem only appears for (p^) but not for (r^), 
so we find the situation inconclusive. 





bb{lS) 


cc{lS) 


bb{lP) 


bb{2S) 


V 


0.26 


0.43 


0.25 


0.24 


yF)(GeV-i) 


1.2 


2.2 


2.1 


2.9 



TABLE I: Estimates for v = \J {p^) /m? and \J (r^) for the Heavy Quarkonium states. For the 
bb{2S) state the number we give for v is quite uncertain. 



The results of this and the following section have been obtained solving the Schroedinger 
equation numerically. We have performed a series of tests of the numerical solutions. As 
we have already mentioned, the case A^ = with z/^ = oo corresponds to the Coulomb case. 
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We have checked the numerical solution against the known analytical result in this case. 
For a general and we have also computed (p^) either directly (in momentum space) 
or through the equality {{E — V{r))) = {p'^/m). Finally, we have also checked the wave- 
function at the origin, either by direct computation (taking the smallest point at which 
the wave function has been computed and checking for stability) or through the equality 
|0nKO)P =W(4vr)(V'(r))„; (see Ref. Q)- 

IV. Ml TRANSITIONS 

In this section we compute the Ml radiative transitions for the low lying bottomonium 
and charmonium states. 

A. T(15) ^ r?fe(15)7 

Our central value for Tr(is)^r]b(is)'r is obtained using Eq. f|TT]) with N = 3, u = 1.5 GeV, 
and i^f = }A- = 0-7 GeV. For we take the values of the T(IS') and ^^^(IS') masses from 



the PDG 



0^ 



In table [TTl we show the size of the different contributions to Tr{is)-^r]i,{is)'y 
The 0{as) and 0{a^) corrections are evaluated at the mass scale. The 0{as{m)) corrections 
in the RS' and on-shell scheme are equal. Renormalon effects first appear at 0{a^{m)) 
and make this expansion more convergent. Yet, as we have taken an small value of Uf, the 
0{a'^{m)) term is still large. There are not 0{v) corrections. The (9(t'^) correction can be 
evaluated at different orders in A^, and for different values of the factorization scale. One 
can easily deduce its size by multiplying Fig. Ob by -5/3 times the LO result. The value 
quoted in Table HTl for the 0{v'^) correction has been obtained for = 3 and u = 1.5 GeV. 
An almost identical value is obtained if one takes u to be the scale of minimal sensitivity. 
Actually, one also obtains a quite similar value if one takes the scale of minimal sensitivity of 
the N = 3, Ur = computation. The great advantage of using z/^ = 0.7 GeV versus z/^ = oo 
is that the u dependence becomes very mild, so it is not a source of uncertainty and one can 



We note though that there is a recent determination of the r]]y{lS) mass which is around 10 MeV lower 



46[ than the PDG value. If such value is confirmed our numbers should be reescaled accordingly (as they 



enter to the third power the effect is important), which can be trivially done. 
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give sensible predictions for the central values. We note that, depending on the order N, 
minimal sensitivity scales may not show up, as we can see in Fig. [3] for other values of 
and/or z/^. Therefore in some cases such prescription may not give a meaningful result and 
the series still be convergent. 





LO 


0(«s) 


0{al) 


0{v^) 




V X 0(^2) 


<5r (eV) 


14.87 


1.29 


0.73 


-1.71 


0.15 


-0.45 



TABLE II: The leading and subleading contributions to rT{i5)-s-»yt,(i5)7- ^c-^^ ^'^^ numbers are 
some error estimates. 

If we sum all the contributions of Table |TT] we obtain 15.18 eV, which is quite close 
to the LO 14.87 eV value. This is due to the strong cancellation between the and v 
corrections. The main source of uncertainty comes from higher order terms. Due to the 
strong cancellation between the as and v terms, we feel that to add a power of v to the overall 
correction would underestimate the error. Instead, we take the v x 0{v'^) contribution in 
Tab. [TTlas our estimate of the subleading correction, as it is the biggest possible contribution. 
Such term alone produces an error of order 3%. This error is much bigger than the error one 
would obtain only considering scale variations (see Fig. [7]), or if we do the evaluation with 
N = 2 instead of = 3 (see again Fig. [Tj), or than the error associated to variations of Uf. 
All these errors are associated to higher order effects. We do not include those in order to 
avoid double counting. The only other source of theoretical error that we include is the one 
due to Nm (for the evaluation of this error we also take into account the correlation with 
the bottom mass value). Beside the theoretical error we also include the error associated to 
the QCD parameters, even though its size is typically smaller than the theoretical error. For 
as we take the variation as{Mz) = 0.118 ± 0.001 [7|. For the variation of the MS bottom 
mass we take mi,{mb) = 4.19 ± 0.03 GeV. In summary, we obtain the following result for the 
different error contributions: 

rriis)^Mis), = 15.18 ±0.45(O(t;3)):°:J2(A^„J;0:04(c,J;0;20(^_) eV, (29) 
which after combining the errors in quadrature reads 

rT(i5H,,,(i5)7 = 15.18(51) eV. (30) 
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0(a) 
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O(a^) 



FIG. 7: Plot ofT-x(is)^r}i,{iS)'y for the case of the bottomonium ground state using the static potential 
V^) at different orders in perturbation theory: N = 0,1,2,3 with Ur = = 0.7 GeV. The 
horizontal line is our central value and the yellow band our final error estimate. 



This corresponds to a branching fraction of 2.9 x 10 Eq. fl3UI) is bigger than result obtained 



in Ref. 



18| (~ {k^/mf X 2.5 keV, see Ref. |47|) but compatible within errors 



B. J/^l>{lS) ^ 7?e(lS)7 

Our central value for T ji^(is)^nc{is)'^ is obtained using Eq. (ITT]) with N = ?>, u = l.b GeV, 
and I'f = Yx = 0.7 GeV. For we take the values of the J/ip{lS) and ric{lS) masses from 
the PDG 7|. In table llllt we show the size of the different contributions to Tj/^(^is)^Ti^(is)-y 
The 0{as) and Ola"^) corrections are evaluated at the mass scale. The 0{as{m)) corrections 
in the RS' and on-shell scheme are equal. Renormalon effects first appear at 0{a'^{m)) 
and make this expansion more convergent. Yet, as we have taken an small value of Uf, the 
0{ag{m)) term is still large. There are not 0{v) corrections. The correction can be 

evaluated at different orders in A^, and for different values of the factorization scale. One 
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can easily deduce its size by multiplying Fig. Hlb by -5/3 times the LO result. The value 
quoted in Table UTTl for the correction has been obtained for = 3 and v = 1.5 GeV. 

Unlike in Sec. IIV Al in this case there are not scales of minimal sensitivity. The use of a 
finite Ur significantly diminishes the factorization scale dependence of the result. Yet we also 
observe that a large scale dependence remains for small scales. Therefore, the value we take 
and quoted in Table UTTl for the correction corresponds to = 3 and u = 1.5 GeV, as 

we feel that smaller values of u may yield unrealistic results. 





LO 


0(as) 


0{al) 


0{v^) 


as X 0{ag) 


V X 0{v^) 


6r (keV) 


2.34 


0.33 


0.16 


-0.71 


0.05 


-0.30 



TABLE III: The leading and subleading contributions to ^ j/^{is)^T]c{iS)'y 
are some error estimates. 



The last two numbers 



If we sum all the contributions of Table IIIII we obtain 2.12 keV, which is quite close 
to the LO 2.34 keV value. This is due to the strong cancellation between the as and v 
corrections. The main source of uncertainty comes from higher order terms. Due to the 
strong cancellation between the as and v terms, we feel that to add a power of v to the overall 
correction would underestimate the error. Instead, we take the v x 0{v'^) contribution in 
Table [1111 as our estimate of the subleading correction, as it is the biggest possible individual 
term. This term alone produces an error of order 15%. This error is much bigger than 
the error one would obtain only considering scale variations (see Fig. |8]), or if we do the 
evaluation with N = 2 instead of A^ = 3 ((see again Fig. |8]). It is also bigger than the error 
associated to variations of Uf. All these errors are associated to higher order effects. We do 
not include those in order to avoid double counting. The only other source of theoretical 
error that we include is the one due to Nm (for the evaluation of this error we also take into 
account the correlation with the charm mass value). Beside the theoretical error we also 
include the error associated to the QCD parameters, even though its size is typically smaller 
than the theoretical error. For as we take the variation as{Mz) = 0.118 ± 0.001 j?]. For the 
variation of the MS charm mass we take mdmc) = 1.25 ± 0.04 GeV (see for instance j4l|). 
In summary, we obtain the following result for the different error contributions: 



rj/^(i5)^,.(i5h = 2.12±0.30(O(i;3))l[J:^^(Ar„);[;:[;^(as);[i:i^(mMs) keV, 



(31) 
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FIG. 8: Plot 0/ rj/^(-x5)_j.^^(i5)^ using the static potential V-^} at different orders in perturbation 
theory: N = 0,1,2,3 with = i^f = 0.7 GeV. The horizontal line is our central value and the 
yellow hand our final error estimate. 



which after combining the errors in quadrature reads 



J/t/.{15)^7y,{15)7 



2.12(40) keV. 



(32) 



This corresponds to a branching fraction of 2.28 x 10 ^. 
We can now compare with previous determinations of this decay. As in Ref. 



48j, we 



summarize the conaparison in Fig. [9l Unhke in this reference we do not include the values 
obtained in Refs. [3, [2^ assigned to potential models. They correspond to the LO com- 
putation in our notation (see Table [In|) . The difference with our value is (mainly) due to 
the different value of the charm mass. Note that in our case the charm mass is not a free 
parameter but fixed by the value of the MS mass. We could still vary the RS' mass by 
changing z//, but this effect would be compensated by the 0(as, v) effects. We now compare 
with the EFT computation of Ref. 18|. It is equivalent to ours setting N = Vj. = 00, 
and eliminating the 0{a1) corrections. When we do so we can get agreement with their 
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FIG. 9: Comparison of different theoretical and experimental predictions for V jj^^^^ 



number if we, as they do, set a very low value for the factorization scale (there are minor 
differences coming from the values of the heavy quark masses used). We see a very strong 
scale dependence in this regime. In this paper we restrict to values of v where we get stable 
results after the resummation of the soft logarithms. This produces much bigger numbers, 
which however get reduced by increasing N . Either way, it should be emphasized that both 
results are perfectly compatible within errors. We also refer to Fig. [H] for comparison with 
the existing experimental numbers 49N51|. and other theoretical predictions using either dis- 



persion relations/sum rules 52|-l54j or lattice simulations 48|, |55|. Our result is compatible 
with most of them within errors. We can discriminate very low values of the decay and start 
to have tensions with the Crystall Ball determination. 



C. P-wave decays 

We now compute the P-wave decays for n = 2 bottomonium (though they could end up 
being of academic interest because of the very small energy differences). In this case we have 
several decays (see Eq. f lT^ ). The difference among them are spin factors, which differently 
weight the fp' matrix element (there are also important differences for k-^ depending on the 
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decay). From the physical point of view the situation is similar to the two previous sections, 
as the wave function still has a single node, though more weighted at somewhat smaller 
scales. Our central value for the decays Th^(_ip)_^^^^(^ip)^, r?,^(ip)^XM(i^')7 and T^^^(^ip)^ht{iP)-y 
are obtained using Eq. ( 11^ with N = 3, u = 1.5 GeV, and Uf = Ur = 0.7 GeV. They 
are compatible with the numbers obtained in [isl if we account for the different and a 
trivial missprint (keV — )■ eV). For k.y we take the masses of the different P-wave states from 
the PDG |7|. In table HVl we show the size of the different contributions. The 0{as) and 

corrections are evaluated at the mass scale. The 0{as{m)) corrections are equal in 
the RS' and on-shell scheme. Renormalon effects first appear at 0{a'^{m)) and make this 
expansion more convergent yet, as we have taken an small value of z//, the (9(as(m)) term is 
still relatively large. There are not 0{v) corrections. The 0{v'^) correction can be evaluated 
at different orders in N, and for different values of the factorization scale. One can easily 
deduce its size by multiplying Fig. Ob by -5/3 times the LO result. Unlike in Sec. IIV Al 
in this case there are not scales of minimal sensitivity. The use of a finite Ur significantly 
diminishes the factorization scale dependence of the result. Yet, we observe that an strong 
scale dependence remains for small scales. Therefore, the value we quote in Table HVl for the 

correction corresponds to = 3 and u = 1.5 GeV, as we feel that smaller values of 
u may yield unrealistic results. Note that, unlike the 0{as) corrections, this contribution is 
differently weighted for each decay. Therefore, properly weighted differences of these decays 
may yield absolute determinations of the (p^)2i matrix element. In any case, the relative sign 
between the and 0{as{m)) corrections produces cancellations between these terms. 

The magnitude of this cancellation depends on the specific decay mode, but they are large 





LO 


0(as) 




0{v^) 


as X ©(a^) 


V X 0{v^) 


5r/i6(ip)^xf,o{ip)7 (^^) 


0.895 


0.078 


0.044 


-0.054 


0.009 


-0.013 


6Th,(iP)^^^^(iP)^ (eV X 10-3) 


8.86 


0.77 


0.43 


-1.08 


0.09 


-0.27 




0.113 


0.010 


0.006 


-0.011 


0.001 


-0.003 



TABLE IV: The leading and subleading contributions to Tfi^(^ip^^^^^(^ip'j^, r/jj^(ip)_j.^j_^(ip)^ and 
Xb2{^P)^h,,{iP)'Y respectively. The last two numbers are some error estimates. 
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FIG. 10: Plot o/r;,^(ip)^^^Q(ip)^, ^hi(iP)^XbiiiPh ^'^'^ Xb2iiP)^ht(iPh using the static potential 
V^} at different orders in perturbation theory: N = 0,1,2,3 with v^. = = 0.7 GeV. The 
horizontal line is our central value and the yellow band our final error estimate. 

in all cases. 

In order to estimate the errors we proceed analogously to the two previous sections. We 
take the v x (9(t>^) contribution of Table HVl as our estimate of the subleading correction, as 
it is the biggest possible individual term. The only other source of theoretical error that we 
include is the one due to (for the evaluation of this error we also take into account the 
correlation with the bottom mass value). Beside the theoretical error we also include the 



error associated to the QCD parameters, even though its size is typical 



theoretical error. For ag we take the variation as{Mz) = 0.118 ± 0.001 I?]. For the variation 
of the MS bottom mass we take mb{mh) = 4.19 ± 0.03 GeV. In summary, we obtain the 



y smaller than the 
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following result for the different error contributions for the three decays: 

r.,(iP)^;,,o(iP), = 0.962 ± O.O13(0(t;3))+°;°29(Ar„);O;OO}(a3);°:°}^(mHs) eV , (33) 
r.,(iP)^;,,(iP), = 8.99 ± 0.27(O(t;3))+0;^6(iV„);0:0^(a,);°:}2(mMs) x IQ-^ eV , (34) 
rx..(iP)^h.{iP)7 = 0.118 ± 0.003(O(t;=^))+°:°°^(iV^);°:°°°(«s);nS^(mMs) eV , (35) 

which after combining the errors in quadrature read 

r,,(iP)^^,„(iP), = 0.962(35) eV, (36) 
r;,,(iP)^^,^(iP)^ = 8.99(55) X 10-3 eV, (37) 
rx..{iPH/..{iP)7 = 0.118(6) eV. (38) 

The errors are heavily dominated by theory. They are much bigger than the error one would 
obtain only considering scale variations (see Fig. [TOj) . or if we do the evaluation with N = 2 
instead of = 3 ((see again Fig. [TOj) . They are also bigger than the error associated to 
variations of Uf. 



D. T(25) ^ r/;,(25)7 

We now compute the T(2S') rjh{2S)'~^ decay. We use Eq. (fTTj) with n = 2, which depends 
on (p^)2o- We observed in Fig. |6]that this object was not convergent in A^. Therefore the 
results of this section should be taken with some caution. 





LO 


0(as) 


0{al) 


0{v'^) 


5T (eV) 


0.640 


0.056 


0.031 


-0.059 



TABLE V: The leading and subleading contributions to rx(25)^77j(25)7- 



Our central value will be obtained using Eq. fllip with N = 3, u = 1.5 GeV, and 
Uf = Ur = 0.7 GeV. For we take the values of the T(25') mass from the PDG [7]. For 
the mass of the rif,{2S) we take the recent value obtained by Belle 46|] for definiteness. 
Nevertheless, we should remark that there is a different value obtained using CLEO data 
561] (if so our numbers can be trivially reescaled accordingly). In table |Vl we show the size of 
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the different contributions to Tr{2S)-¥ri,,{2S)-y The 0{as) and O^a"^) correction are evaluated 
at the mass scale. The 0{as{m)) corrections are equal in the RS' and on-shell scheme. 
Renormalon corrections first appear at 0{a'^{m)) and make this expansion more convergent 
yet, as we have taken an small value of Uf, the 0{a'^{m)) term is still large. There are not 
0{v) corrections. The 0{v'^) correction can be evaluated at different orders in A^, and for 
different values of the factorization scale. One can easily deduce its size by multiplying Fig. 
Ob by -5/3 times the LO result. The value quoted in Table IVl for the correction has 

been obtained for = 3 and u = 1.5 GeV. For the 2S bottomonium state the use of a 
finite Ur, in particular Ur = 0.7, does not significantly improve over the Ur = oo computation. 
The factorization scale dependence is still significant and the convergence bad. Therefore, 
conservatively we take the ^) term as our estimate of the error associated to higher order 
corrections. This term alone produces an error of order 10%. For the rest of the errors we 
proceed as in the previous sections. Overall, we obtain 

rT(25)^..(25h = 0.668 ± Om9{0{v'))t'ofotiNrnno^^^^^^ eV , (39) 

which after combining the errors in quadrature reads 

rT(25H^,(25)7 = 0.668(60) eV. (40) 

In Fig. [11] we compare this result with the scale variation of the evaluation of the decay for 
different values of A^. Note that our error is much bigger than the one would assign from 
the factorization scale dependence, or from the difference between different evaluations. 
We believe an error analysis only based on those would underestimate the error. 

E. 2S ^IS-f decays 

The experimental situation of the 25—7-15' radiative transitions has improved signif- 
icantly over the last years. Whereas for the 2^5o — )■ 1^5i7 decay there is still no data 
available, this is not so for the 2^5i — ?■ 1^5o7 decay, for which the PDG [7] quotes the value 
[3.9 ±1.5] X 10~^ for the decay branching fraction. This translates into the following value 
for the decay 

r?gH,,(i5h = 12.5(49) eV. (41) 
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FIG. 11: Plot of rx(2S)-i.r;j,(2S)7 for the case of the bottomonium ground state using the static 
potential V^) at different orders in perturbation theory: N = 0,1,2,3 with Ur = Uf = 0.7 GeV. 
The horizontal line is our central value and the yellow band our final error estimate. 

This number comes from js^] BaBar (branching fraction [3.9±l.l{stat)^Ql{syst)] x 10~^) and 
updates the previous upper bound branching fraction < 0.5 x 10^^ ( or rx(2S)^%(i5')7 < 0.016) 
produced by CLEOIII 

On the theoretical side the 25 — ?■ IS* radiative transitions are different from the previous 
transitions considered before. Now we only know the leading nonvanishing order, see Eqs. 
( IT2|) and ( 1131) . which scales as ~ {k^/m'^)v^. It depends on the expectation values of p^, 

and Vs2{r) among different states {n = 1 and n = 2), which we have not studied so far 
(note that for those matrix elements we can only fix their relative sign but not the absolute 
one). Moreover the V52(r)/m^ potential is modulated by the Wilson coefficient D^g2 which 
resums the large logarithms associated to the heavy quark mass. 

The warning qualifications that we made in the previous section may also apply here, as 
the decay depends on the dynamics of the 2S bound state, for which we have observed prob- 
lems of convergence in for (p^)2o- Nevertheless it is worth repeating that the observables 
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FIG. 12: Plot o/ rx(25)^r?5(i5)7 using the static potential at different orders in perturbation 

theory: = 0, 1, 2, 3 with Vf = 0.7 and = oo. The dashed line corresponds to no resumming the 
hard logarithms: Dg2 g = as(z^). The blue band corresponds to the experimental value. 

we are sensitive to now are different, and, therefore, worth exploring them. In Fig. [T^ we 
show the theoretical predictions for Tr{2S)^'nb(is)-y after approximating the static potential 
at different orders in working at z/^ = oo and Uf = 0.7 GeV (see solid lines). In other 
words, we just add an r-independent constant to the static potential. We actually see a 
nicely convergent pattern for the decay, the magnitude of which decreases quite significantly 
as we increase (by around an order of magnitude) and approaches the experimental value. 

In order to understand this result it is convenient to study the magnitude of the different 
terms that contribute to Eq. f|T2|) . We display the terms inside the brackets of Eq. f|T2|) in 
Fig. [13] with z/j. = oo and Uf = 0.7 GeV (see dashed lines). Note that they are O^v"^) ~ 0.06 
up to prefactors. We observe a very nice convergence pattern for the io('"^)2o associated 
term. The convergence of the io(j5^)2o term is less good, and even less for io(Ks2(^))2o (for 
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V (GeV) 

FIG. 13: Plot of 1st, 2nd and 3rd term inside the brackets of Eq. using the static potential 

V^) at different orders in perturbation theory: N = 0,1,2,3. The dashed lines correspond to the 
Vf = 0.7 GeV and Vr = oo case, and the continuous lines to the Vf = = 0.7 GeV case. 

scales below 2.5 GeV). In any case, there is a very strong cancellation between the different 
terms in the decay. This makes the total sum of these terms smaller than the magnitude 
of each of them. The bulk of this effect is independent of the factorization scale and gets 
strongly magnified as we increase (see again the solid lines of Fig. [T2l) . Therefore it does 
not seems to be a numerical accident for an specific or factorization scale. 

If we switch off the resummation of the hard logarithms and work at = with z/^ = oo, 
our computation is equivalent to the analysis performed in Ref. jisl. In this reference 
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a very large value for the decay was obtained. We show our equivalent computation as 
the dashed line in Fig. [121 If we set z/ = 1 GeV, the value used in this reference, we 
obtain rT(2S)-i>r)i,(is)7 — 0.6 keV (~ 0.649 keV if we use the mass mf, = 4730 MeV used in 
that reference). Therefore, the introduction of the hard logarithms are crucial to make the 
decay transition width smaller for = at small scales. As we increase N this effect is 
less important and the decay width gets small irrespectively of the resummation of hard 
logarithms (yet the final value may change by a factor 3 or 4, specially at small scales). 
At this stage we would like to emphasize that the computation of the decay shows a nicely 
convergent pattern in A^, as we can see in Fig. [12l rapidly approaching the experimental 
number. 

As in previous sections we can try to improve the previous results by exactly incorporating 
the correct asymptotic short distance behavior of the static potential in the solution of the 
Schroedinger equation. Typically the convergence is accelerated and the factorization scale 
dependence greatly diminishes. We show the behavior of the different contributions to the 
decay in Fig. [13] (see solid lines). There is a strong cancellation between the second and 
third term in Eq. ( |T2|) (compare the solid lines of Fig. [13lb and Fig. [I3lc), whereas the first 
term is almost constant (see the solid lines of Fig. [T3la). We show the result for the decay in 
Fig. [Hlwith i>r = 0.7 GeV, where we also compare with the z/,. = oo case, and experiment. 
Note how this figure corresponds to a zoom of Fig. [12], as the scale dependence is much 
smaller, as well as the size of the corrections. For A^ = we still have some scale dependence 
but it basically vanishes for A^ > and u > 2 GeV. Actually, the results are very stable to 
scale variations (with a nice plateau for u > 2 GeV) and to the value of A^. In order to get 
these results the resummation of the hard logs play a crucial role, specially a low u. We also 
emphasize that the final numbers compare quite favourable with experiment. This is by far 
nontrivial, as there has been more than one order of magnitude reduction with respect the 
original numbers obtained with a pure Coulomb potential. 

Therefore we dare to give a value for, and assign errors to, the decay width. In order to 
produce our final numbers we proceed analogously to the previous sections. In Table [12] we 
give the prefactor, the different matrix elements for A^ = 3, = 1.5 GeV and z/,. = z/j = 0.7 
GeV, as well as the total decay width. In order to estimate the error associated to subleading 
effects in f , we multiply the biggest of the three f ^ contributions by v, instead of multiplying 
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0(a) 



O(a^) 



O(a^) 



FIG. 14: Plot of Tx{2S)^T]t,iiS)j using the static potential at different orders in perturbation 

theory: N = 0,1,2,3 with Uf = = 0.7 GeV (solid lines) and Vf = 0.7 GeV, v,. = oo (dashed 
lines). The blue hand corresponds to the experimental value and the yellow hand to our theoretical 
prediction. 

by V the sum of the three contributions, as we can not guarantee that the cancellation 
between different terms takes place at higher orders. The structure of the error estimate 
would then be 



"J- T(2S)^%(1S)7 



B[{Av' + 5v^f - {Av'f] ~ B{Av')2 x 5v' - 0.005 keV, (42) 



where A is an small number and 5v^ ~ t> x 0{v'^). We check the reliability of this error 
estimate by replacing the theoretical masses that appear in the third term in Eq. ( IT2l) by 
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Prefactor (keV) 




A{f) 


A{Vs2) 


r(eV) 


10.3342 


0.022 


0.039 


-0.042 


6.3 



TABLE VI: The prefactor, the terms inside the brackets of Eq. U^) . and the total decay width 
rT{25)^%(15)7- 



the physical ones (as our result is very sensitive to this term). This effect is subleading in v 
and produces a shift with respect to the central value of order STr(2S)-^r]t(is)'r — 0.006 keV. 
We take this number (which is quite similar to the number obtained in Eq. fH^ ) as our 
estimate of the higher order uncertanties. We do not dwell further in the analysis of the 
higher order uncertanties, as our main error will come from an strong dependence in A^m^- 
We also compute the error associated to as and m^. Summarizing all the errors we obtain 



.(th) 



0.006 ± 0.006(0(t;^))™(iV^);HlJU«s);noo("^Ms) keV. (43) 



If we combine all the errors in quadrature our final number reads 



.(th) 

T(2S)->%(1S)7 



6+2^ eV 

"-06 ^ ^ ■ 



(44) 



The error is completely dominated by theory. It completely covers the experimental predic- 
tion. Note that the same error is obtained using the scale variation of the = 3, z/^ = oo 
result, which does not depend on A^^ but is typically less precise. 

Overall we conclude that Tr{2S)-*r]b{is)'r is relatively suppressed with respect its natural 
size by a very large cancellation between the (p^) and {Vs^) term. This makes the total matrix 
element smaller^. The fact that it enters as magnifies this effect. A confirmation of this 
picture would come from the evaluation (and experimental determination) of r^j(2S)^T(i5')7, 
which we expect to be much larger because the relative sign between these two terms changes. 
Actually this is what we find, as one can see in Fig. [151 On the other hand, for this decay, 
there is no convergent pattern in A^. Therefore, we do not dare to make any error analysis, 
and only estimate the decay to be around r|^*'|25)^T(i5)7 ~ 



^ We could reduce the dependence on Nm by increasing vf (and Vr)- The price one would pay is a stronger 

dependence on v. 
^ And an ideal place to measure 1 10(^^)20 1- 
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We can compare Eq. f H^ with the recent lattice simulation of Ref. [59|]. As our com- 
putation is 0{v^) we should compare with their 0{v^) result. In matrix element units it 
corresponds to their number 0.080(5) in Table II (the experimental number is 0.035(7)). Our 
central value corresponds to 0.025!l;;o;g25- Nevertheless, a proper comparison would require 
the incorporation of the renormalization group improved Wilson coefficient, Dg2 in the 
lattice analysis^. If we switch it off in our analysis our result gets strongly scale dependent 
and we can get agreement with their results for an scale of around 1 GeV. At this respect 
we can not avoid mention that we expect some dependence on the lattice spacing of the 
NRQCD matrix elements as, in general, it is not possible to obtain the continuum limit 
for them. In any case, we now face an interesting situation: In Ref. jsQl agreement with 
experiment was only obtained after the inclusion of the 0{v^) operators (again using tree 
level Wilson coefficients). Note that this implies a complete breakdown of the v expansion 
for bottomonium, as the 0{v^) correction would be as important as the 0{v'^) term. On the 
other hand our picture is different and is possible to obtain agreement with experiment with 
an Olv"^) computation (and the help of the renormalization group at small scales). 



V. CONCLUSIONS 



We have computed the magnetic dipole transitions between low lying Heavy Quarkonium 
states in a model independent way. We have used the weak-coupling version of pNRQCD 
with the static potential exactly incorporated in the LO Hamiltonian. The precision we have 
reached is k^/m^ x (9(ag,f^) and k^/m^ x (9(t>^) for the allowed and forbidden transitions 
respectively. Large logarithms associated to the heavy quark mass scale have also been re- 
summed. The effect of the new power counting was found to be large and the exact treatment 
of the soft logarithms of the static potential made the factorization scale dependence much 
smaller. The convergence for the hh ground state was quite good, and also quite reasonable 
for the cc ground state and the hh IP state. For all of them we gave solid predictions, which 

^ This has not been done so far. We stress that this effect could be quite important. Actually one only has to 
incorporate the very same logs that we are incorporating here, as the leading logs are scheme independent. 
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O(a^) 

O(a^) 
0(a) 

FIG. 15: Plot o/ r^j^(25)^x{i5)7 using the static potential V^gV at different orders in perturbation 
theory: N = 0,1,2,3 with Ur = Vf = 0.7 GeV. 

we summarize here: 



rT(lS)^r,b(lS)7 


= 15.18(51) eV, 


(45) 


rj/i/.{15)^»7c(lS)7 


= 2.12(40) keV, 


(46) 


^ hi{lP)^XboO-P)l 


= 0.962(35) eV, 


(47) 




= 8.99(55) X 10-^ eV, 


(48) 




= 0.118(6) eV. 


(49) 



For the 25* decays the situation is less conclusive. The 0{v'^) correction of the T(2S) — )■ 
r]b{2S) 7 decay suffered of bad convergence in producing relatively large errors for our 
prediction (see Eq. fHOj) ). Some of the Oiy"^) matrix elements of the ?7b(25') — > T(IS') 7 decay 
also suffered of this bad convergence. This made impossible to give a reliable error estimate 
for this transition, as such terms correspond to the leading (and only known so far) order 
expression (moreover they should be squared in the decay). The situation is completely 
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different for the T{2S) — )■ rib{lS)'y transition. The reason is that the problematic 
matrix elements appear in a different combination for this decay, so that they cancel to a 
large extent. This led to a nicely convergent sequence in A^, where the resummation of the 
hard logarithms played an important role. Our final figure was 

T(2S)^%(1S)7 ~ "-06 ^ V • 

This number is perfectly consistent with existing data, so that the previous disagreement 
with experiment for the T{2S) — t- rib{lS) 7 decay fades away. 

The error of the above figures is dominated by theory, in most cases by the lack of 
knowledge of higher order effects. The determination of the origin and nature of those 
effects may significantly diminishes the errors. Typically, they may come from loop effects, 
so it may happen that they effectively count as 0{asv'^), implying smaller errors. In any 
case, let us note that the static potential becomes steeper as we increase N. Therefore, the 
transfer energy between the heavy quarks is bigger and the effective alpha and radius of 
the bound state becomes smaller than what one would deduce from a pure LO Coulomb 
evaluation. This means that the weak-coupling approximation works better than one would 
expect a priori for those systems. These are good news for weak-coupling analysis of the 
properties of the lowest-lying heavy quarkonium resonances. 

We have not incorporated the error associated to in our final numbers. Therefore, 
strictly speaking, our figures are theoretical predictions ofT/k^. In some cases the associated 
error would be small. Yet, we have chosen to work this way since, for some decays, the 
experimental value of k^ is still uncertain. It is trivial for the reader to introduce such error 
if willing to do so. 

We have also profited to compute some expectation values like the electromagnetic radius, 
(r^), or (p^). We find (r^) to be nicely convergent in all cases, whereas the convergence of 
(p^) is typically worse. We have found that (p^) is more or less constant with n, and 



a/ (r^) is more or less linear with n. This is the same behavior one finds with a logarithmic 
potential. Since the early days of heavy quarkonium it is well known from potential models 
that such potential effectively describes the spectrum of the bottomonium and charmonium 



system 



44( 1 . We find rewarding that the QCD potential can simulate such behavior after the 



inclusion of the running of as 
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The computation of (r^) and (p^) (and the binding energies) also yield a very nice con- 
firmation of the renormalon dominance picture. This predicts in the on-shell scheme that, 
on the one hand, the binding energy should diverge with but, on the other, (r^) and (p^) 
should produce convergent sequences in A^. We have observed this effect in full glory in our 
analysis. 
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